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Abstract. A weak-lensing analysis of a wide- field F-band 
image centered on the cluster Abell 1942 has uncovered 
the presence of a mass concentration projected ~ 7 ar- 
cminutes South of the cluster center. From an additional 
wide-field image, taken with a different camera in the I- 
band, the presence of this mass concentration is confirmed. 
A statistical analysis, using the aperture mass technique, 
shows that the probability of finding such a mass concen- 
tration from a random alignment of background galaxies 
is 10~^ and 4 x 10"** for the V- and /-band image, respec- 
tively. No obvious strong concentration of bright galaxies 
is seen at the position of the mass concentration, but a 
slight galaxy number overdensity is present about 1' away 
from its center. Archival ROSAT-HRI data show the pres- 
ence of a weak extended X-ray source near to the mass 
concentration, but also displaced by about 1' from its 
center, and very close to the center of the slight galaxy 
number concentration. 

From the spatial dependence of the tangential align- 
ment around the center of the mass concentration, a rough 
mass estimate can be obtained which depends strongly 
on the assumed redshift of the lens and the redshift 
distribution of the background galaxies. A lower bound 
on the mass inside a sphere of radius 0.5h~^ Mpc is 
1 X 10^'*/i~^AfQ, considerably higher than crude mass esti- 
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mates based on the X-ray data; shifting the lens to higher 
redshift increases both the lensing and X-ray mass esti- 
mates, but does not resolve the mass discrepancy. 

Concerning the nature of the mass concentration, no 
firm conclusion can be obtained from the available data. If 
it were a high-redshift cluster, the weak X-ray flux would 
indicate that it had an untypically low X-ray luminosity 
for its mass; if the X-ray emission were physically unre- 
lated to the mass concentration, e.g., coming from a rel- 
atively low-redshift group which shows up in the number 
density of galaxies, this conclusion would be even stronger. 

Since the search for massive halos by weak lensing en- 
ables us for the first time to select halos based on their 
mass properties only, it is possible that new types of ob- 
jects can be detected, e.g., halos with very little X-ray 
and/or optical luminosity, should they exist. The mass 
concentration in the field of A1942 may be the first ex- 
ample of such a halo. Possibilities to establish the nature 
of this mass concentration with future observations are 
briefly discussed. 

Key words: Cosmology: dark matter, gravitational lenses 



1. Introduction 

The abundance of clusters of galaxies as a function of mass 
and redshift provides one of the most sensitive cosmolog- 
ical tests (e.g., Richstonc et al. 1992; Bartelmann et al. 
1993). In particular, in a high-density Universe, the abun- 
dance of massive clusters strongly decreases with redshift, 
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so that the existence of a few massive high-redshift clus- 
ters can in principle rule out an f2o = 1 model (e.g., Eke 
et al. 1996; Bahcall & Fan 1998). 

The reliability of the test depends on the detection ef- 
ficiency and selection effects in existing samples of clusters 
whose understanding may be critical. Currently, clusters 
are selected either by their optical appearance as overden- 
sities of galaxies projected onto the sky and/or in color- 
magnitude diagrams, or by their X-ray emission. Both se- 
lection techniques may bias the resulting sample towards 
high-luminosity objects, i.e. they would under-represent 
clusters with high mass-to- (optical or X-ray) light ratio. 
Furthermore, the observed properties have to be related 
to their mass in order to compare the observed abundance 
to cosmological predictions. The usual procedures consist 
in assuming a dynamical and/or hydrostatic equilibrium 
state as well as the geometry of the mass distribution, 
which in general may be questionable and fairly poorly 
justified from a theoretical point of view. 

Indeed, whereas cosmological theories have made great 
progress in their ability to predict the distribution of dark 
matter in the Universe, either analytically or numerically 
(e.g., Lacey & Cole 1993; Jenkins et al. 1998), the lu- 
minous properties of matter are much more difficult to 
model. For example, to relate the X-ray data of a cluster 
to its mass, a redshift-dependent luminosity- temperature 
relation needs to be employed (sec Borgani et al. 1999 
and references therein), in the absence of a detailed un- 
derstanding of the physics in the intra-cluster gas. It would 
therefore be of considerable interest to be able to define a 
sample of 'clusters' - or more precisely, dark matter ha- 
los - which can be directly compared with the predictions 
coming from N-body simulations. 

Weak gravitational lensing offers an attractive possi- 
bility to detect dark matter halos by their mass properties 
only. A mass concentration produces a tidal gravitational 
field which distorts the light bundles from background 
sources. Owing to their assumed random intrinsic orienta- 
tion, this tidal field can be detected statistically as a co- 
herent tangential alignment of galaxy images around the 
mass concentration. A method to quantify this tangential 
alignment was originally introduced by Kaiser et al. (1994) 
to obtain lower boimds on cluster masses, and later gener- 
alized and proposed as a tool for the search of dark mat- 
ter halos (Schneider 1996). This so-called aperture mass 
method can be applied to blank field imaging surveys to 
detect peaks in the projected density field. Combining halo 
abundance predictions from Press & Schechter (1974) the- 
ory with the universal density profile found in N-body sim- 
ulations (Navarro et al. 1997), Kruse & Schneider (1999) 
estimated the number density of dark matter halos de- 
tectable with this method (with a signal-to-noise thresh- 
old of 5) to be of order 10 deg~^, for a number density of 
30 galaxies/ arcmin^, and depending on the cosmological 
model. These predictions were confirmed (Reblinsky et al. 



1999) in ray-tracing simulations (Jain et al. 1999) through 

numerically-generated cosmic density fields. 

In this paper, we report the first detection of a dark 
matter halo not obviously associated with light, using the 
above-mentioned weak lensing technique. Using a 14' x 14' 
deep V-hand image, obtained with MOCAM at CFHT, we 
aimed to investigate the projected mass profile of the clus- 
ter Abell 1942 on which the image is centered. We foimd a 
highly significant peak in the reconstructed mass map, in 
addition to that corresponding to the cluster itself. This 
second peak, located about 7' South of the cluster cen- 
ter, shows up in the alignment statistics of background 
galaxy images with a significance > 99.99%, as obtained 
from Monte-Carlo simulations which randomized the ori- 
entation of these background galaxies. An additional deep 
7-band image, taken with the UII8K at CFHT, confirms 
the presence of the mass peak. No obvious large overdcn- 
sity of galaxies is seen at this location, implying either a 
mass concentration with low light-to-mass ratio, or a halo 
at substantially higher redshift than A1942 itself. Finally, 
an analysis of an archival ROSAT/HRI image of A1942 
shows, in addition to the emission from the cluster, a 3.2- 
a detection of a source with position close to the peak 
in the projected mass maps; though this weak detection 
would be of no significance by itself, the positional coinci- 
dence with the 'dark' clump suggests that it corresponds 
to the same halo, and that it may be due to a high-redshift 
{z k, 0.5) cluster. 

The outline of the paper is as follows: in Sect. 2 we 
describe the observations and data reduction techniques, 
as well as the measurement of galaxy ellipticities which 
we employed. The aperture mass statistics is briefly de- 
scribed in Sect. 3.1 and applied to the optical data sets, 
together with a determination of the peak detection sig- 
nificance. Properties of the mass concentration as derived 
from the optical data sets and the X-ray data are dis- 
cussed in Sects. 3.2 and 3.3, respectively, and a discussion 
of our findings is provided in Sect. 4. We shall concentrate 
in this paper mainly on the 'dark' clump; an analysis of 
the mass profile of the cluster A1942 and the reliability 
of mass reconstruction will be published elsewhere (van 
Waerbeke et al., in preparation) 

2. Summary of optical observations and image pro- 
cessing 

The V- and /-band observations were obtained at the 
prime focus of CFHT with the MOCAM and the UH8K 
cameras, respectively. Both observing procedures were 
similar, with elementary exposure time of 1800 seconds 
each in V and 1200 seconds in /. A small shift of 10 arc- 
seconds between pointings was applied in order to remove 
cosmic rays and to prepare a super-flatfield. 

The V"-band images were obtained during an observing 
run in dark time of June 1995 with the 4K x AK mosaic 
camera MOCAM (Cuillandre et al 1997). Each individual 
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Fig. 1. Tlie geometry of tlie optical data used in this paper. The left-hand side shows the area of the K-band MOCAM field 
(square) and the 7- band UH8K-chip3 data (rectangle). The framed regions are 3'.3 x 3'.3 cutouts around the cluster center of 
A1942 and around our 'dark clump' candidate. These regions are zoomed in on the right-hand side. The 'dark clump' region is 
centered around q( J2000) = 14'' 38™ 22.59"; 5(J2000)=03° 32' 32.22". 



chip is a X 2K LORAL CCD, with 0'.'206 per pixel, so 
the total field-of-view is 14' x 14'. Nine images have been 
re-centered and co-added, to produce a final frame with a 
total exposure time of 4h30min. The seeing of the coadded 
image is 0'.'74. 

The /-band images were obtained with the x 
mosaic camera UH8K (Luppino, Bredthauer & Geary, 
1994). Each individual chip is a 2if x AK LORAL CCD, 
also with 0'.'206 per pixel, giving a field-of-view of 28' x 28'. 
The final centered coadded image resulting from 9 sub- 
images has a total exposure time of 3h and a seeing of 
0'.'67. The V- and /-band images have been processed in 
a similar manner, using standard IRAF procedures and 
some more_specific ones developed at CFHT and at the 
TERAPIX0 data center for large-field CCD cameras. None 
of these procedures had innovative algorithms, so there is 
basically no difference in the pre-processing and process- 
ing of the MOCAM and UH8K images. For the present pa- 
per, we use only Chip 3 of the UH8K /-band image which 

^ http://terapix.iap.fr 



contains the cluster A1942, and the additional mass con- 
centration discussed further below. Fig. 1 shows the CCD 
images from both fields and their relative geometry. 

A first object detection and the photometry have 
been performed with SExtractor2.0.17 (Bertin & Arnouts 
1996). The MOCAM field has been calibrated using the 
photometric standard stars of the Landolt field SAllO 
(Landolt 1992), and the UH8K field was calibrated using 
the Landolt fields SA104 and SAllO. The completeness 
limit is y = 26 and / = 24.5. 

The lensing analysis was done with the imcat software, 
based on the method for analysing weak shear data by 
Kaiser, Squires & Broadhurst (1995), with modifications 
described in Luppino & Kaiser (1997) and Hoekstra et al. 
(1998; hereafter HFKS98). This method is based on cal- 
culations of weighted moments of the light distribution. 
Imcat is specifically designed for the measurement of el- 
lipticities of faint and small galaxy images, and their cor- 
rection for the smearing of images by a PSF, and for any 
anisotropy of the PSF which could mimic a shear signal. 
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Fig. 2. For all objects detected with the imcat method from the MOCAM frame (left panel) and the UH8K chip3 (right panel), 
the magnitude is plotted as a function of half-light radius rh, measured in pixels. Objects containing saturated pixels have been 
removed from the plots. The magnitudes are in an arbitrary system. In both cases we can clearly identify a prominent sequence 
of stellar objects at about rh = 2.2 for MOCAM, and r-h = 1.75 for UH8K-chip3. 



These corrections are employed by the relation 

X = X° + ^^7 + ^'^>, (1) 

where x is the observed image ellipticity (defined as in, 
e.g., Schneider & Seitz 1995), x" is the ellipticity of the 
unlensed source smeared by the isotropic part of the PSF, 
is the response tensor of the image ellipticity to a shear, 
and P^™ is the response tensor to an anisotropic part of 
the PSF, characterized by p. These tensors are calculated 
for each galaxy image individually. Since the expectation 
value of x" in (1) is zero, one obtains an unbiased estimate 
of the shear through 

7={P'')-'[X-P'"'P] ■ (2) 

(7 is in reality an estimate for the reduced shear 7/(1 — k) 
which reduces to the shear if k <^ 1.) The PSF anisotropy 
in our images is fairly small and regular over the field. We 
selected bright, unsaturated stars from a size vs. magni- 
tude plot (see Fig. 2) and determined their ellipticities. As 
Fig. 3 shows, the stellar ellipticity changes very smoothly 
over the fields so that its behaviour can be easily fit with a 



second-order polynomial (see also Fig. 4). With these poly- 
nomials we performed the anisotropy correction in (1). We 
follow the prescription of HFKS98 for the calculation of 
P'', and used the full tensors, not just their trace-part, in 
(2). 

The current version of imcat does not give information 
about the quality of objects; for this we produced a SEx- 
tractor (version 2.0.20) catalog containing all objects that 
had at least six connected pixels with 1-cr above the local 
sky background. From this catalog we sorted out all ob- 
jects with potential problems for shape estimation (like be- 
ing deblended with another object or having a close neigh- 
bour). This included all objects with FLAGS> 2 (internal 
SExtractor flag) . The remaining catalog was matched with 
the corresponding imcat catalog, using a maximum posi- 
tional difference of three pixels, and keeping only those 
objects for which the detection signal-to-noise of imcat 
was > 7. 

This procedure left us with 4190 objects {V > 22.0) for 
the MOCAM and 1708 objects (/ > 21.0) for the 7-band 
chip3. With these final catalogs all subsequent analysis 
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Fig. 3. The cllipticity fields for stars for the V-band MOCAM field (left panel) and the /-band UH8K-chip3 containing the 
cluster. Both fields show a smooth variation and can be easily modelled by a low-order polynomial. The maximal ellipticity is 
about 5% for the MOCAM and 8% for the UH8K. 



was done. We note that we did not cross-correlate the 
MOCAM and UH8K catalogs; hence, the galaxies taken 
from both catalogs will be different even in the region of 
overlap. Due to the different waveband used for object se- 
lection, the redshift distribution of the background galax- 
ies selected on the MOCAM and the UH8K-chip3 frame 
can be different. 

3. Analysis of the 'dark' clump 

3.1. Weak tensing analysis 

From the image ellipticities of 'background' galaxies, we 
have first reconstructed the two-dimensional mass map 
of the cluster field from the MOCAM data, using the 
maximum-likelihood method described in Bartelmann et 
al. (1996) and independently, the method described in 
Seitz & Schneider (1998). The resulting mass maps are 
very similar, and we show the former of these only. 

In the left panel of Fig. 5, we show the resulting mass 
map with the (mass-sheet degeneracy) transformation pa- 
rameter A chosen such that (k) = (see Schneider & Seitz 



1995) , together with contours of the smoothed number 
density of bright galaxies. In general, this number density 
correlates quite well with the reconstructed surface mass 
density. As can be seen, a prominent mass peak shows up 
centered right on the brightest cluster galaxy. 

In addition to this mass peak, several other peaks are 
present in the mass map. Such peaks may partly be due to 
noise coming from the intrinsic image ellipticities and, to 
a lesser degree, to errors in the determination of image el- 
lipticities. In order to test the statistical significance of the 
mass peaks, we used the aperture mass method (Schneider 

1996) . 

Let U (?9) be a filter function which vanishes for -d > 9, 
and which has zero mean, Jq di) ^U{'&) = 0. Then we 
define the aperture mass Map(i?) at position ■& as 

Map(^?) = / d^t?' k(i? + U{\^'\) . (3) 

J\'d'\<e 

Hence, Map(i?) is a filtered version of the density field k; 
it is invariant with respect to adding a homogeneous mass 
sheet or a linear density field, and is positive if centered 
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Fig. 4. The left panels show the raw imcat ellipticitics from bright, unsaturated foreground stars in our fields (upper pan- 
els: MOCAM field; lower panels: UH8K chip). The right panels show the ellipticities after they have been corrected with a 
second-order polynomial as described in the text. The rms of the ellipticities after correction is typically 0.015. 



on a mass peak with size eomparable to the filter scale 9. 
The nice feature about this aperture mass is that it can 
be expressed directly in terms of the shear, as 

MapW= / , dV7t(i?';i?)g(|i?'|) (4) 

Jl-d \<e 

(Kaiser et al. 1994; Schneider 1996), where the filter func- 
tion Q{i}) = 2i9"2 J^^ 1?' U{i9') - t/(i9) is determined 
in terms of U{'&), and vanishes for •& > 9. The tangential 
shear 7t(i?'; i?) at relative position with respect to i? is 
defined as 

7t(^?'; 1?) = -7ee[7(i? + 1?') e'^"^'] , (5) 

where ip' is the polar angle of the vector . In the case 
of weak Icnsing (k ^ 1), the observed image ellipticities 7 
from (2) are an unbiased estimator of the local shear, and 
so the aperture mass can be obtained by summing over 
image ellipticities as 

Mi^{^) =—j2iuWQm-^\), (6) 



where the sum extends over all N galaxy images with 
positions di which are located within 9 of ■d, and the tan- 
gential component 7tj('«?) of the image ellipticity relative 
to the position is defined in analogy to 71. In general, 
M^p(i?) is not an unbiased estimator of Map(i?) since the 
expectation value of 7 is the reduced shear, not the shear 
itself. However, unless the aperture includes a strong mass 
clump where k is not small compared to unity, M^p will 
approximate Map closely. But even if the weak-lensing ap- 
proximation breaks down for part of the aperture, one can 
consider the quantity M^p(i?) in its own right, represent- 
ing the tangential alignment of galaxy images with respect 
to the point i9. This interpretation also remains valid if the 
aperture is centered on a position which is less than 9 away 
from the boundary of the data field, so that part of the 
aperture is located outside the data field, in which case 
M^p(i?) will not be a reliable estimator of Map(i?). 

In order to determine the significance of the peaks in 
the mass map shown in Fig. 5, we have calculated M!^p 
on a grid of points t? over the data field, for four values 
of the filter scale 9. Then, we have randomized the posi- 
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Fig. 5. The figure shows mass reconstructions and galgixy number density from the MOCAM field (left panel) and the 
UH8K-chip3 (right panel). The white contours show k = 0.03, 0.05, 0.07, 0.1, 0.12, 0.15, 0.17 and 0.2. For the reconstruc- 
tion the shear was smoothed with a Gaussian of cr = 40" width. The black contours show the smoothed galaxy distribution 
from all galaxies brighter than V = 21.0 and / = 20.0 (the smoothing kernel here was a Gaussian with cr = 20"). 



tion angles of all galaxy images, and calculated M^^ on 

the same grid for these randomized realizations. This has 
been repeated iVrand times. Finally, at each grid point the 
fraction v of randomizations where M^p is larger than the 
measured value from the actual data has been obtained; 
this fraction (which we shall call 'error level' in the follow- 
ing) is the probability of finding a value of M^p at that 
gridpoint for randomly oriented galaxy images, but with 
the same positions and ellipticities as the observed galax- 
ies. 

Fig. 6 displays the contours of constant p, for different 
filter radii, varying from 80" to 200". As can be seen, the 
cluster center shows up prominently in the I'-map on all 
scales. In addition, two highly significant peaks show up, 
one at the upper right corner, the other ~ 7' South of 
the cluster center, close to the edge of the MOCAM field. 
We have verified the robustness of this Southern peak by 
using SExtractor ellipticities instead of those from imcat, 
and found both the cluster components and the Southern 
peak also with that catalog (although it should be much 
less suited for weak lensing techniques). 

After these findings, we obtained the UH8K /-band 
image, on which both the cluster and the Southern mass 
peak are located on Chip 3. The mass reconstruction from 
galaxy images on Chip 3 are shown in the right panel of 
Fig. 5, from which we see that the cluster and this South- 
ern mass peak also show up. Repeating the aperture mass 



statistics for Chip 3, we obtain the error levels as shown 

in Fig. 7; again, this Southern peak shows up at very high 
significance. Whereas the third peak in the significance 
maps (considering the two larger filter scales) from Chip 
3, about halfway between cluster and the Southern compo- 
nent and slightly to the West, is also quite significant and 
is also seen in the corresponding MOCAM map (and most 
likely also corresponds to a mass peak, though a highly 
elongated one for which the aperture mass is less sensi- 
tive), we shall concentrate on the Southern peak, which 
we call, for lack of a better name, the 'dark clump'. 

In fact, as can be seen from Figures 1 and 5, this mass 
peak does not seem to be associated with any concen- 
tration of brighter galaxies. This could mean two things: 
either, the mass concentration is in fact associated with 
little light, or is at much higher redshift than A1942 itself. 

Concentrating on the location of the dark clump, we 
determined the probability distribution po{M!^p) for the 
value of M^p, obtained from 2 x 10® randomizations of 
the galaxy orientations within 160" of the dark clump. 
This probability distribution is shown as the solid (from 
MOCAM) and dashed (from Chip 3) curve on the left 
of Fig. 8. These two distributions are very well approxi- 
mated by a Gaussian, as expected from the central limit 
theorem. The value of Mj^p at the dark clump is 0.0395 
for MOCAM, and 0.0283 for Chip 3. The fact that these 
two values are different is not problematic, since for Chip 
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Fig. 6. The four panels show the significance v (sec text) of tlic Af^p maps of the MOCAM field. Wc chose A'^rand = 5000, 
the black contours mark areas with v = \, 10, 30/5000 and the white contours v = 100, 180, 260/5000. The filter scales are 80" 
(upper left panel), 120" (upper right panel), 160" (lower left panel) and 200" (lower right panel). For the larger scales the cluster 
components and the dark clump are detected with a very high significance 



3, the whole aperture fits inside the data field, whereas it 
is partially outside for MOCAM; hence, the two values of 
M^p measure a different tangential alignment. Also, since 
the two data sets use galaxies selected in a different wave- 
band, their redshift distribution can be different, yielding 
difTcrcnt values of the resulting lens strength. The proba- 
bility that a randomization of image orientations yields a 
value of Mj^p larger than the observed one is ~ 10~^ for 
the MOCAM field, and 4.2 x 10"^ for Chip 3. 

Next we investigate whether the highly significant 
value of M^p at the dark clump comes from a few galaxy 
images only. For this, the sample of galaxy images inside 
the aperture was bootstrap rcsamplcd, to obtain the prob- 
ability pboot(-M^^p) that this resampling yields a particular 
value of M^p. This probability is also shown in Fig. 8. The 
probability that the bootstrapped value of M' is nega- 



tive is 3.8 X 10-"^ for Chip 3, and < IQ-^ for the MOCAM 
peak. 

The radial dependence of the tangential image elliptic- 

ity is considered next. Fig. 9 shows the mean tangential 
image ellipticity in annuli of width 20", both for the MO- 
CAM and the UH8K data centered on the dark clump. 
The error bars show the 80% probability interval obtained 
again from bootstrapping. It is reassuring that the radial 
behaviour of (7t) is very similar on the two data sets. 
In fact, owing to the different wavebands of the two data 
fields and the fact that the aperture does not fit inside the 
MOCAM field, this agreement is better than one might 
expect. The mean tangential ellipticity is positive over a 
large angular range; except for one of the inner bins (for 
which the error bar is fairly large), (7t) is positive in all 
bins for Q ^ 150". This figure thus shows that the large 
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Fig. 7. The same as Fig. 6 for the UH8K-chip3 7-band data. The filter scales are, from left to right: 80", 120", 160" and 200". 
A1942 and the dark clump are also detected here with a very high significance 



and significant value of M^p at the dark clump is not domi- 
nated by galaxy images at a particular angular separation. 

3.2. Properties of the dark clump 

We now investigate some physical properties of our dark 

clump candidate. Wc first argue that it is very unlikely 
for our object to lie at a redshift higher than 1. For our 
magnitude limit of 24.5 in the I band we expect approx- 
imately 30 galaxies/ (1')^. We used approximately half of 
them (see Sec. 2) as putative background galaxies for our 
analysis. The median of simulated redshift distributions 
that extend the CFRS data (Lilly et al. 1995) to fainter 
magnitude limits (Baugh, Cole & Frenk 1996) is at about 
0.7 — 0.8. If we assume that all our galaxies lie in 
the extreme tail of these distributions, then z = 1.0 rep- 
resents a good upper limit for the redshift of our clump. 
However, the lensing analysis of the high-redshift cluster 
MS1054— 03 (Luppino & Kaiser 1997) may provide an in- 
dication for a somewhat larger mean source redshift. 

Next we use Fig. 9 to obtain a crude estimate of the 
mass of this object. Although the tangential shear appears 
to be fairly small close to the center position of the clump, 
there is a region between ~ 50" and ~ 150" where the tan- 
gential shear is clearly positive and decreases smoothly 
with radius. If we describe the mass profile by an isother- 
mal sphere, its velocity dispersion cr„ would be given by 



rameters, this becomes 



(7) 



where the product 7t 9 would be independent of 9 for an 
isothermal sphere model, and the final term is the ratio 
lens-source to observer-source distance, averaged over the 
background galaxy population. Introducing fiducial pa- 



= 1135 



Tioo 



0.06 V 3 (£»ds/£'s 



■ km/s , 



(8) 



where 7100 is the tangential shear 100" from the mass 
center. Alternatively, we can express this result in terms of 
the mass within a sphere of radius R, M{< R) = 2alR/G; 
for example, within R = 0.5h~^ Mpc, we find 



M(< 0.5/1-^ Mpc) = 2.9 X lO^"* h-^MQ 



1 



0.06 3{Di,/D,) ■ 
(9) 

Whereas this model is quite crude, the largest uncertainty 
in quantitative mass estimates comes from the unknown 
redshift of the dark clump and the unknown redshift dis- 
tribution of the background galaxy population. The mass 
is a monotonically increasing function of the lens redshift, 
and depends very strongly on the assumed mean source 
redshift, in particular for values of ^ 0.5. 

With the / band data wc now estimate the light coming 
from the dark clump. For this we created a SExtractor cat- 
alog counting every connected area with at least 3 pixels 
0.5-(T above the sky backgroimd as a potential object. The 
flux of all these objects (except from obvious stars) in a 
circle of 100" radius around the clump center was summed 
up. Wc did the same in 32 control circles around 'empty' 
regions in the other UH8K chips. It turned out that the 
flux within the clump region is compatible with the mean 
flux of the control annuli, i.e., there is no ovcrdcnsity of 
light at the position of the dark clump. So we took the 1-a 
fluctuation of the fluxes in the control circles as a reason- 
able upper limit for the light coming from the dark clump. 
For converting the flux into a total / band magnitude we 
assumed that we are dominated by elliptical galaxies, us- 
ing K corrections for this galaxy type calculated with the 
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Fig. 8. Probability distributions for M^p, with the aperture centered on the peak position of the dark clump. Solid (dashed) 
curves correspond to the MOCAM (Chip 3) data set. For an aperture of 160", the left of the two curves shows the probability 
distribution po{M^p) for values of M^p obtained by randomizing the position angles of the galaxy images. These two curves 
nearly coincide. The two curves on the right-hand side show the probability distribution Pboot(Afap) obtained from bootstrap 
resampling of the galaxy images inside the aperture. The two vertical lines show the measured values of M^p 



latest version of the Bruzual & Chariot stellar population 
synthesis models for the spectrophotometric evolution of 
galaxies (Bruzual fc Chariot 1993). Prom the total / band 
magnitude we derived a bolomctric magnitude and a bolo- 
metric luminosity using standard approximations. With a 
lower limit for the mass and an upper limit for the lumi- 
nosity wc can give lower limits for the mass to light ratio 
of our object. This is shown in Fig. 10 for different source 
redshift distributions and two cosmologies. We see that 
the EdS imivcrsc gives fairly high M/L estimates in com- 
parison to a n = 0.3, A = 0.7 model. When we assume a 
redshift of 2 « 0.8 for our clump we obtain a lower limit 
of M/L « 300 in the A cosmology. This is a conservative 
lower limit which could be lowered significantly only if one 



assumes that the redshift distribution of the faint galaxies 
extends to substantially higher redshift. 

As the dark clump has a mass characteristic of mas- 
sive clusters it is of interest to search for X-ray emission 
associated with it. 

3.3. The X-Ray data analysis 

A1942 was observed by the ROSAT HRI in August 1995. 
The total integration time was 44,515 s. We retrieved the 
X-ray images from the public archive and reduced them 
using ESAS, Snowden's code especially developed for the 
analysis of extended sources in ROSAT data (Snowden et 
al 1994; Snowden & Kuntz 1998). 

The region showing a significant peak in the weak lens- 
ing reconstructed mass map is within the field of view of 
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Fig. 9. Mean tangential image ellipticity in independent bins of width 20" around tlie dark clump, triangles show the mean, 
solid (dashed) error bars the 80% error interval obtained from bootstrapping, using the MOCAM (Chip 3) data. For better 
display, the points and error bars are slightly shifted in the 6 direction. 



the HRI image of A1942. Wc have searched for X-ray emis- 
sion in this area. First of aU, wc have refined the astrom- 
etry in the X-ray image matching X-ray point sources to 
objects in om deep optical images. The astrometric offset 
from the original instrument coordinates is 3.5" . There is a 
significant X-ray emission peak centered at 14"^ 38™ 22.8**, 
3° 33' 11" (J2000.0). This position is 60" away from the 
weak lensing mass peak. The X-ray source is detected at 
the 3.2-cr level using an aperture of 30" radius. Although 
the number of counts detected is low, its distribution is in- 
consistent with a point-like source, showing a profile elon- 
gated along the NW-SE direction that is broader than the 
instrumental PSF. 

Wc have measured the source count-rate using con- 
centric circular apertures centered on the X-ray emission 
peak. We obtain a count-rate of 7.4±2.5 x 10~S~^ within 
a circular aperture of 45" radius. The counts still increase 



somewhat at larger radii but the measurement is much 
noisier given the uncertainty in the sky determination. 
The total flux is thus approximately 10-30% larger than 
the above value. We convert the count-rate into a flux as- 
suming an incident spectrum of T = 3 keV and a local 
hydrogen column density of Nh = 2.61 x 10^^ cm~^. The 
resulting unabsorbed flux is 3.4± 1.2 x 10"-'^^ erg cm~^ s""'^ 
in the 0.1-2.4 keV band. We have also fitted a standard 
beta profile (Cavaliere & Fusco-Fermiano 1978) to the az- 
imuthaUy averaged radial profile. We obtain best values 
for the core radius and beta parameter (slope decline at 
large radii) of 15" and 0.80, respectively, although these 
values are quite uncertain given the low total number of 
counts. 

The X-ray luminosity depends on the redshift of the 
source. Assuming an incident spectrum at the detector of 
T = 3 keV [T = 3(1 -|-2:)keV at the source], the rest-frame 
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Fig. 10. Estimate of the lensing mass (upper left panel), an upper bound for the luminosity of the lens (upper right panel), and 
a lower limit on the mass-to- light ratio (lower panel), as a function of assumed lens redshift. All estimates are for an aperture 
size of 100". The solid, short dashed and long dashed curves show the MjL ratio in an EdS universe for (zs) = 0.8, (2a) = 0.9 
and (zs) = 1.0. The dotted, dot-short dashed and dot-long dashed curves show the same in an = 0.3 ,A = 0.7 universe. Wc 
have assumed a redshift distribution oc exp {-{zlzaf^] for the source galaxies; hence {zs) ~ l.Szo. A value of 7100 = 0.06 
was assumed. 



X-ray luminosity in the 0.1-2.4 kcV band would range from 
1.9 ± 2.5 X 10''2/i-2 erg s^^ if the redshift is the same as 
that of A1942 (z = 0.223) to 3.5 ± 0.5 x lO^^/i-^ erg s'^ 
\iz = 1.0 {qo = 0.5). 

We have also made a crude estimate of the mass of 
the system. On the one hand, if we assume an X-ray 
luminosity-temperature relation (e.g., Reichart et al 1999, 
Arnaud & Evrard 1999) and a temperature- mass relation 
(e.g., Mohr et al 1999), we can get mass estimates at a 
0.5/1-1 ]y|-pp i-adius from 1.5 x IO^^/i-^Mq at 2; = 0.223 to 
1.6 X IO^^H^^Mq at z = 1 ((7o = 0.5). We can also assume 
a beta profile, fixing the core radius and the beta param- 
eter, and compute the normalization necessary to obtain 
the observed flux at the measured radius. Then we can 



integrate the profile to obtain the gas mass. If we further 
assume a gas fraction, wc can also obtain a total mass 
estimate. If we take the values obtained from our previ- 
oTis fit of the X-ray surface brightness profile, we get total 
masses at a radius of 0.5h-^ Mpe, of 9.2 x W^'^h^^MQ 
at z = 0.223 and 2.3 x IO'^^H-^Mq at z = 1 {qo = 0.5). 
Note the difference of a factor of 1.5 and 7 compared to 
the previous estimates. This gives an indication of the er- 
rors involved. If instead we were to use typical values of 
the core radius and beta parameter of most clusters of 
galaxies (e.g., Tc = 0.125/1^^ Mpc and (3 = 2/3) the mass 
estimates would be approximately a factor 3 larger and 
closer to the estimates using standard correlations. 
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Fig. 11. For the UH8K-chip3 field, surface mass density (black) and X-ray (white) contours arc plotted. The surface density 
contours are the same as in Fig. 5, whereas the X-ray contours correspond to 1.5 x 10~^, 1.6 x 10^'', 2.0 x 10^'', 3.0 x 10~^ and 
4.0 X 10~® counts/s/pixel. The cluster A1942 itself is clearly seen in X-rays, centered on the brightest cluster galaxy. In addition, 
extended X-ray emission near the dark clump is detected 



Although we have presented quantitative values for the 
mass of the system based on the X-ray data, these should 
be taken only as informative given the assumptions and 
errors involved. Our main point in presenting these esti- 
mates is to show that this system has the X-ray properties 
of a galaxy group if it is at the same redshift of A1942. 
The lensing shear signal measured would then be too large 
for such a group unless it had a remarkable unusually high 
mass-to-X-ray light ratio. It seems more plausible that the 
system is a more massive cluster of galaxies at a higher 
redshift if the X-ray and lensing signal do indeed come 
from the same source, although the X-ray derived mass 
is still lower than the one obtained from the shear signal. 
The small angular scale X-ray core radius (larger physical 
scale if at larger redshift) and the lack of bright galaxies 
also point towards the same conclusion. 



As an alternative, the X-ray emission may be unrelated 
to the dark clump, but associated with the small galaxy 
number overdensity projected near it, as scon from the 
black contours in the right-hand panel of Fig. 5. In that 
case, both the local enhancement of the galaxy density 
and the X-ray emission may be compatible with a group 
of galaxies, rather than a massive cluster, as indicated by 
the weak lensing analysis. 

4. Discussion and conclusions 

Using weak lensing analysis on a deep high-quality wide- 
field V-band image centered on the cluster Abell 1942, we 
have detected a mass concentration some 7' South of the 
cluster. This detection was confirmed by a deep 7-band 
image. No clear overdensity of bright galaxies spatially 
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associated with this mass concentration is seen; therefore, 
wc termed it the 'dark chimp'. A shght ovcrdcnsity of 
galaxies is seen ~ 1' away from the mass center of the dark 
clump, but it is unclear at present whether it is physically 
associated with the mass concentration. Archival X-ray 
data allowed us to detect a 3.2-a X-ray source near the 
dark clump, separated by 60 arcseconds from its peak; 
it appears to be extended. The X-ray source is spatially 
coincident with the slight galaxy overdensity. 

We have estimated the significance of the detection of 
this mass peak, using several methods. For the l/-band im- 
age, the probability that this mass peak is caused by ran- 
dom noise of the intrinsic galaxy ellipticities is ~ 10~^; a 
similar estimate from the I-band image yields a probability 
of ~ 4 X 10^*. Thus, the mass peak is detected with ex- 
tremely high statistical significance. A bootstrapping anal- 
ysis has shown that the tangential image alignment is not 
dominated by a few galaxy images, as also confirmed by 
the smooth dependence of the tangential shear on the an- 
gular separation from its center. Whereas these statistical 
tests cannot exclude any systematic effect during obser- 
vations, data reduction, and ellipticity determination, the 
fact that this dark clump is seen in two independent im- 
ages, taken in different filters and with different cameras, 
make such systematics as the cause for the strong align- 
ment highly unlikely. Although we have accounted for the 
slight anisotropy of the PSF, the uncorrected image ellip- 
ticities yield approximately the same result. 

A simple mass estimate of the dark clump shows it to 
be truly massive, with the exact value depending strongly 
on its redshift and the redshift distribution of the faint 
background galaxies. The mass inside a sphere of ra- 
dius 0.5h~^ Mpc is ^ 10^'^h~^MQ, if an isothermal sphere 
model is assumed; if the lens redshift is larger, this lower 
mass limit increases, by about a factor 2 for z ^ 0.5 and 
a factor of about 10 for z 1. In any case, this mass esti- 
mate appears to be incompatible with the X-ray flux if the 
dark clump corresponds to a 'normal cluster', at any red- 
shift. We therefore conclude that the mass concentration, 
though of a mass that is characteristic of a massive cluster, 
is not a typical cluster. This conclusion is independent of 
whether the X-ray emission is physically associated with 
the dark clump or not. 

The lack of an obvious concentration of galaxies near 
the mass peak has been transformed into an upper limits 
of the luminosity associated with the mass concentration, 
and therefore into a lower limit of the mass-to-light ratio. 
This M/ L limit depends again strongly on the redshift dis- 
tribution of the faint galaxies, as well as on the assumed 
clump redshift. Whereas values for M/L as low as ^ 200 
(in solar units) are theoretically possible if the clump has 
a redshift in excess of unity, the corresponding mass be- 
comes excessively and unrealistically large; for more rea- 
sonable redshifts ^ 0.8, M/L ^ 450 for an Einstein-de 
Sitter Universe, and M/ L ^ 300 for a low density flat Uni- 
verse. We would like to point out, though, that estimates 



of the M/ L-ratio quoted in the literature practically never 

assume a A-dominated cosmology, so that the M/L ratio 
quoted above for the low-density Universe cannot be di- 
rectly compared to literature values. 

We can only speculate about the nature of this dark 
clump. As argued above, a normal cluster seems to be 
ruled out, owing to the lack of bright X-ray emission. 
Whereas the estimated X-ray luminosity can be increased 
by shifting the putative cluster to higher redshifts, the 
corresponding lens mass also increases with Zd, in a way 
which depends on the redshift distribution of the source 
galaxies. The spatial coincidence of the slight galaxy over- 
density and the X-ray emission, both ~ 1' away from the 
mass center of the dark clump, may best be interpreted as 
a galaxy group or weak cluster at relatively low redshift 
and not associated with the dark clump. 

The dark clump itself may then be a mass concentra- 
tion with either low baryon density or low temperature, 
or both. For example, it may correspond to a cluster in 
the process of formation where the gas has not yet been 
heated to the virial temperature so that the X-ray lumi- 
nosity is much lower than expected for a relaxed cluster. 
The fact that the tangential shear decreases towards the 
center of the mass clump may indeed be an indication of 
a non-relaxed halo. 

Further observations may elucidate the nature of this 
mass concentration. Deep infrared images of this region 
will allow us to check whether an overdensity of IR- 
selected galaxies can be detected, as would be expected 
for a high-redshift cluster, together with an early-type 
sequence in the color-magnitude diagram. A deep image 
with the Hubble Space Telescope would yield a higher- 
resolution mass map of the dark clump, owing to the large 
number density of galaxies for which a shape can be mea- 
sured, and thus determine its radial profile with better ac- 
curacy. Images in additional (optical and IR) wavebands 
can be used to estimate photometric redshifts for the back- 
ground galaxies. In conjunction with an HST image, one 
might obtain 'tomographic' information, i.e., measuring 
the lens strength as a function of background source red- 
shift; this would then yield an estimate of the lens redshift. 
The upcoming X-ray missions will be considerably more 
sensitive than the ROSAT HRI and will therefore be able 
to study the nature of the X-ray source in much more de- 
tail. And finally, one could seek a Sunyaev-Zel'dovich sig- 
nature towards the dark clump; its redshift-independence 
may be ideal to verify the nature of a high-redshift mass 
concentration. 

But whatever the interpretation at this point, one must 
bear in mind that weak lensing opens up a new channel 
for the detection of massive halos in the Universe, so that 
one should perhaps not be surprised to flnd a new class 
of objects, or members of a class of objects with unusual 
properties. The potential consequences of the existence of 
such highly underluminous objects may be far reaching: 
if, besides the known optical and X-ray luminous clusters. 
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a population of far less luminous dark matter halos exist, 
the normalization of the power spectrum may need to be 
revised, and the estimate of the mean mass density of the 
Universe from its luminosity density and an average mass- 
to-light ratio may change. We also remind the reader that 
already for one cluster, MS1224, an apparently very high 
mass-to-light ratio has been inferred by two completely 
independent studies (Fahlman et al. 1994; Fischer 1999). 
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